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^ ; Abstract 

I A nonlinear dynamics of the thermal and electromagnetic instabilities of critical state in type II 

' superconductors has been analysed taking into account the effect of dissipation and dispersion. An 

existance of a nonlinear running waves describing a final stage of evolution of the thermomagnetic 
pH I instability in superconductors is demonstrated analitically. 
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5-1 ■ The evolutions of the of the thermal (T) and electromagnetic {E^ H) perturbations is 

^ i' 

^ ' described by the nonlinear equation of the thermal conductivity [1] 

c/3 ; 

■ dT , , ->-^ 

|; z/— = V[«:VT]+jE, (1) 

^ ' by the Maxwel's equations 
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O : rotE = — — , 2 

I— —I, c dt 
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■ rotH = —J (3) 
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^ I and by the equation of the resistive state 

j=Ut,H)+1{E). (4) 

' where z/ = z/(T) and k = n(T) are the heat capacity and the thermal conductivity respec- 

S ' tively; jc is the critical current density and jr is the resistance current density. 

' The above system is essentially nonlinear because the right-hand part of Eq.(l) contains a 

Q I term describing the Joule heat evolution in the region of a resistive phase. Such a set (l)-(4) 

H i of nonlinear parabolic differential equations in partial derivatives has no exact analytical 

[ solution [1]. 

^ ! Let us consider a planar semi-infinite sample (x > 0) placed in an external magnetic field 

<^ ; - dll 

H = {0,0, Hf.) growing at a constant rate = const. According to the Maxwel's equation 

(2), there is a vortex electric field E = {0,Ee,0) in the sample, directed parallel to the 
current density j: E \\ j; where H^, is the amplitude of the external magnetic field and E^ is 
the amplitude of the external electric field. 

A set system decision (l)-(4) in the form of functions from new automodel variable C,{x,t) 
is presented: 
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T^T[C{x,t)], 



(5) 



Substituting (5) in a set system (l)-(4) gives as a result simple differentation the following 
system 



dt 



dT 



(f^dT ^ /dC 
dx^ d^ \dx . 



d^T ^ (dT 



de \d^ 



+ iUT) + jn{E)]E , 



(6) 



d'^^dE U^yd^E_A7i 
dx"^ d^ \dx I d^"^ 



djc dT djn dE 
dT~d^~^dE~d^ 



di 
dt 



(7) 



That a set system (6), (7) at the substation (5) was only function from ^ is required per- 
forming the following conditions: 



(8) 



(9) 



(10) 



where A, C, C are the functions from ^, type which will be determined below. Solving first 
two system equations (8)-(10) we have a relation 



Whence just filows relationship between G and A type 

til 



(12) 



where u is the free constant of integrating an equation (11). From (8) and (9) follows that 
C{x,t) must satisfy single-line equation in private derived 



d^ _ d^ 
dt dx 



single deciding which is function 



^{x,t) ^F{x-ut) 



(13) 



(14) 



Using (14) immediately get, 
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G(0 = 1, A{0^-Fu, C(0 = 0. (15) 

By the transformation the coordinates and time possible to ensure F = 1. Thereby, 
definitively find sought automodel substitution 

^ — X — ut, (16) 

corresponding to solution a type of running wave [2] . 

For the automodelling solution of the form (16), describing a running wave moving at a 
constant velocity v along the x axis, the system of equations (l)-(4) takes the following form 

= (17) 
'^'j, (18) 



c 

E = -H. (19) 

c 

The thermal and electrodynamic boundary conditinos for equations (17)- (19) are as follows: 

dT 

T{C ^ +oo) = To, —iC ^ -oo) = 0, 



E{i ^ +oo) = 0, E{i ^ -oo) = E„ 



(20) 



where Tq is the temperature of the cooling medium. 

Let us consider the Bean- London model of the critical state for the dependence jc{T, H) 
[3] 

je(r)=io[l-a(T-ro)] (21) 

where jo is the equilibtium current density, a is the thermal haet softening coefficient of the 
magnetic flux pinning force. 
A characteristic field dependence oi jr{E) in the region of sufficiently strong electric fields 

2 

(E > Ef) can be aproximated by a piecewise linear function jV ~ (^fE-, where o"/ = ~ 

anHc^/H is the effictive conductivity in the flux flow regime; rj is the viscous coefficient, $o = 
Trhc 

— — is the magnetic flux quantum, f7„ is the conductivity in the normal state; Ef is the 
boundary of the linear area in the volt age- current characteristics of sample. In the flux 
creep regime {E < Ef) the relation between the current density j, and the electric fleld, E, 
is strongly nonlinear [4]. 
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Excluding are variables T(^) and H{^) using Eqs.(17) and (19) , and taking into account 
the boundary conditions (20), we obtain an equation describing the electric field E{^) dis- 
tribution (£^-wave): 



+ 



Attv djr dE Anv^a N{T) - iV(To) 



dE dC 



aE^ 
2k{T) 



= 0, 



(22) 



dE" 



where the dependency T — T ( — j is defined by expression (2), (4) and have the form 



Attv dE' 



(23) 



Here N{T) = ju{T)dT. 



dE. 



The analysis of the phase plane (E, — ) of Eq.(22) shows that there are two equihbrium 
points: Eq ^ 0, T ^ Tq is the stabile node and E ^ E^, T ^ T* ^ T{Ee, 0) is the saddle. 
The sepatetrix, joining these two equilibrium points represents a solution of the shock- wave- 
type with amplitude E^ (see, [7]). The velocity of an E'-wave is determined by the Eq. (22) 
with account a boundary conditions (20): 

^2 



Stt 



N 



To + -[jc{T)+Jrm 



- iV(To) 



(24) 



Using Eqs. (19) and (22) we find the expression for a distribution of the magnetic field H 
in the case of H-wme : 



d^H Attv 
'de^~d^ 



djr 



dE 



dH ^ ^^N{T) - N{To 



E= 



-H 



di 



H- VH 



N{T) 



k{T) 



av 

Yck{T) 



(25) 



(26) 



The condition (9) reflects tha adiabatical character of the wave propagation: the magnetic 
field energy transferred by the wave ensures local heating of the sample in the close vicinity 
of the wave front. Thereby the depending on the external conditions at the sample surfase 
a possible to exist two kinds thermomagnetic waves. 

Note that within the linear voltage-current chatacteristics {E > Ef) in the approximation 
of weak superconductor heating (T — Tq) «Tc — Tq), the heat capacity u and the thermal 
conductivity k depend on temperature profile only slightly. Then in the limiting case r > > 1 
[5] the effect of dispersion on the wave dynamics is negigible. Correspondingly equation (22) 
coincides with the first integral of Burgers equation [6] , which have solution of the form [7] 
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E{z) 



E* 

T 



1 - th^{z - zq) 



(27) 



Here the following dimensionless parameters are introduced E* 



-; where L 



t 



4vrjo 



nf 

2/3 p=-^andz 
Ij 



is the depth of the magnetic field penetration into the superconductor, 



is the time of thermal diffusion, r 



of the thermal Dt — — to the magnetic Dm — 

constant. 



is the parameter describing the ratio 



diffusion coefficient; E,^ — is a 



Using the boundary condition E{z — > — oo) = Eg we can be determined the velocity Ve 



L r E^ 



2tE, 



the front width 



Sz = 16 



rV2 



Ek 

'eI 



1/2 



and the time of a wave propagating inside the sample 

L 



2vf 



(28) 
(29) 
(30) 

10~^ sec for 



sin 

Numerical estimations give ve ^ 1 -r- 10^ — , Sz ~ 10~^ -j- 10~^ and At 

sec 

characteristics values of the physical parameters [5]. 

The system Eqs. (l)-(4) is invariant with respect to an arbitrary translation. Therefore, the 
wave propagation conditions can be found for an arbitrary critical current density dependence 
on T and H. The results can also be obtained for an arbitrary temperature dependence of 
thermophysical parameters u and k of superconducting material and for an arbitrary function 

jr{E). 

The problem of stability of nonlinear shock-wave with respect to small thermal and elec- 
tromagnetic perturbations is studied in [7] . It is also shown that only damped perturbations 
correspond to space-limited solutions, which means that nonlinear wave is stable. 
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